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I. Introduction 



This report describes position estimation procedures that 
are based on models which relate positional uncertainty to 
various measurement and estimation errors. The procedures were 
developed to be used in analyses that relate the effect of posi- 
tional uncertainty on tactical performance through such factors 
as weapon accuracy as well as to be used operationally. 

In the models, positions are on a plane surface (flat earth 
model). Because of this condition, the models are not intended 
for use in situations where the earth's figure is significant. 

In addition, positional errors are determined by independent 
normally distributed random variables with known means, variances 
and covariances. The support for this condition, other than its 
mathematical convenience, is that it has been used by others, 
for example, see Reference 1. To use the models, one is required 
to specify the means, variances and covariances of the error 
random variables. For these models, this can be done by speci- 
fying the systematic errors (biases) and the error magnitudes 
(standard deviations) . 

The procedure that is described in Section II relates position 
estimates that are based on bearings on or from stations to 
station position uncertainty. One application for the model is 
the analysis of the effect of sonobuoy position uncertainty on 
position estimates determined with passive directional sonobuoys . 

The procedure that is described in Section III relates posi- 
tion estimates that are based on lines of position to line of 
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position uncertainty. In addition to estimating the uncertainty 
in position estimates based on celestial observations, the pro- 
cedure could be used to determine error ellipses for LORAN fixes 
if the standard deviation values required by the procedure could 
be obtained. 

The procedure that is described in Section IV can be used 
to combine position estimates from various sources. The proce- 
dure which is based on conditions that should not be too restric- 
tive in most cases provides both a composite position estimate 
and error ellipse. 
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II. station Position Uncertainty and Position Estimates 

The procedure that is described in this section relates 
position estimates that are based on bearings on or from stations 
to station position uncertainty. The procedure is based on a 
model that is an extension of one that is described in Appendix 
1. The model is defined as follows: Each station position error 

is determined by an independent bivariate normal distribution 
with a zero mean vector and a known covariance matrix. Observed 
bearing lines on or from a station are parallel to true bearing 
lines. The distance of each observed bearing line from its true 
bearing line is determined by an independent normal distribution 
with a zero mean and a standard deviation o. 

Because a station's position error is determined by a bivar- 
iate normal distribution, the perpendicular distance s betv/een a 
line at the assumed location of an observed bearing line and the 
observed bearing line is determined by a normal distribution. 

The relation between the bivariate normal distribution that 
describes the station's position uncertainty is indicated in 
Figure 1. In the figure, the positive y-axis direction is north, 
the positive x-axis direction is east, and the origin of the 
coordinate system is at the assumed station position. The x'y'- 
coordinate system is oriented so that the positive x'axis is 
coincident with the major axes of the elliptical contours of 
the bivariate normal distribution that determines the station 
position error and so that the bearing 6 of the positive x'-axis 
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Figure 1. The geometry associated with the determination of Ug , 
the standard deviation of the normal distribution 
that determines the distance between the assumed 
location of an observed bearing line and the observed 
bearing line. The assumed (mean) position of the 
station is at the origin. The ellipse represents a 
contour on the probability density surface of the 
bivariate normal distribution that determines the 
station's position. 
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satisfies the condition; 0° < 6 < 180°. The angle 6 is an 

observed bearing. The x"y "-coordinate system is oriented so 
that the positive y"-axis is in the direction of the observed 
bearing line. As a consequence of these relationships, 

2 2 2 2 2 
o = 0 ^, sin (6-9) + o„, cos (6-9) 

S A y 

2 2 

where 0 , and o , are the elements of the station position 
X y 

error covariance matrix relative to the x ' y ' -coordinate system. 

A procedure for determining position estimates that are 
based on bearings on or from stations at known positions is 
described in Appendix 1 of this report. The procedure is based 
on a model that is equivalent to one that is described in Ref- 
erence 2. The procedure that is described in this section is 
based on a model that is an extension of it. 

The procedure in Appendix 1 of this report is based on a 
model in which a station's bearing error is determined by a 
normal distribution with zero mean (bias) and standard devia- 
tion e. The bearing error is related to the distance on a cir- 
cular arc between a station's true bearing line and observed 
bearing line. The arc is on the circle with its center at the 
station that passes through an initial estimate of an object's 
position. This distance is determined by a normal distribution 
with mean zero and standard deviation o = re where r is the 
range of the initial estimate from the station and the standard 
deviation e is measured in radians. In the model, arc distance 
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is approximated using a first order approximation which in 
effect replaces the circle with its tangent line at the initial 
estimate. As a consequence, the distance u on the tangent line 
between the observed bearing line and the true bearing line is 
determined by a normal random variable with mean zero and stan- 
dard deviation o. The distance u can be expressed in terms of 
w, the distance on the tangent line between the observed bearing 
line and the initial estimate that is also determined by a normal 
random variable with standard deviation a, and v, the distance 
on the tangent line between the true bearing line and the initial 
estimate. And, as shown in Appendix 1 of this report, this 
distance can be expressed in terms of the unknown coordinates 
of the object's position. 

The effect of station position uncertainty is accounted for 
by the distance s between the observed bearing line and the 
assumed observed bearing line. In the model s is determined by 
a normal distribution with mean zero and standard deviation a 

s 

as given above. This approximation is consistent with the first 
order approximation of arc distance. As a consequence of these 
two approximations, all of the bearing lines are replaced by 
lines parallel to the line joining the initial estimate's 
position and the station's assumed position, both of which are 
known positions. The geometry involved in shown in Figure 2. 

The modified relationships resulting from the introduction of 
station position uncertainty are shown in Figure 3. The 
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Figure 2. The geometry of the position estimation model. The 
bearing 3 and range r of the initial estimate from 
the assumed station position have the role of 3 
and r in Appendix 1 of this report • 
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assumed location 
of the observed bearing line 



Figure 3 




The quantities u, v and w correspond to the quantities 
u, V and w in Appendix 1 of Reference 3. The auxiliary 
quantities u' and w' are defined in the figure. The 
replacement of by + Og^ in the procedure in 
Appendix 1 of this report is justified by noting that 
with this substitution u' and w' are equivalent to u 
and w. 



8 



important thing to note is that u' = u - s is determined by a 

normal distribution with mean zero and standard deviation 
2 2 h 

[<? + Og ] that otherwise U' is equivalent to U with 

respect to the procedure in Appendix 1 of this report. As a 
consequence of this, the procedure can be extended to include 
station uncertainty by replacing o by [o + ^ where ever 

it is used. In this case, o = re where r is the range of an 
initial estimate of an object's position from the assumed station 
position and e is the bearing error (standard deviation) in 
radians of the bearings associated with the station. 
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III. Line of Position Uncertainty and Position Estimates 

The procedure that is described in this section relates 
position estimates that are based on lines of position to line 
of position uncertainty. The procedure is based on a model that 
is defined as follows: Lines of position are straight lines. 

Observed lines of position are parallel to true lines of posi- 
tion. The distance of an observed line of position from a true 
line of position is determined by an independent normally dis- 
tributed random variable with known mean and standard deviation. 
Lines of position are specified in terms of a rectangular coor- 
dinate system with the origin at a reference point as shown in 
Figure 4. For celestial navigation, an appropriate choice for 
the reference point would be the assumed position. Since bear- 
ing lines are lines of position, this model differs from the 
model that is described in Section II only in terminology. 
However, operationally the use of the model that is described 
in this section differs in the way the standard deviation of the 
distance of the line from a true line is determined. The stan- 
dard deviation associated with each line of position must be 
specified. If this can be done, the procedure can be used. 

As an example, suppose that the values are for the first 

line of position and the second where > O 2 and 

that the lines of position intersect at a 90° angle. In this 
example, the minimum area confidence (probability) region is 
an ellipse that is centered on the estimated position. And the 
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position 



o 



1st line of position 



initial estimate 
and final estimate 




2nd line of positic 



Figure 4. The geometry when only two lines of position are used 
with the procedure in Appendix 1 of this report. The lines of 
position correspond to bearing lines for observed bearings and 
0^ with respect to the procedure. The procedure determines the 
size and the orientation of elliptical confidence regions of con- 
fidence p, the location of the initial estimate and the location 
of the final estimate. The locations of the estimates correspond 
when only two lines are used. If the lines were obtained from 
sextant observations, the assumed position should be the reference 
position in which case the lines of position would be determined 
by the azimuth anges and and the distances and P 2 . 
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estimated position is at the intersection of the bearing lines. 

2 

For a confidence (probability) of containment of 1 - exp(-k /2) , 
the major axis of the ellipse is coincident with the first bear- 
ing line and it is of length 2k0j^, the minor axis of the ellipse 

is coincident with the second bearing line and it is of length 

2 

2 k 02 and the area of the ellipse is Trk • 
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IV. A Composite Position Estimate 

The procedure that is described in this section is for 
combining position estimates for an object from independent 
sources. It is based on the following model: The rectangular 

coordinates of each position estimate are determined by an inde- 
pendent bivariate normal distribution whose covariance matrix 
is known but whose mean vector is not known. The components of 
the mean vector for each of the distributions are the unknown 
coordinates x and y of the object. This model implies that the 
natural logarithm of the likelihood function for a set of n 
estimates can be expressed as follows: 

n I -1 

log L = K - 1/2 Z (x. - x) T. . (x- - x) 

i=l 

where K is a constant, x^^ is an estimate vector with components x^^ 
and y^, x is the common mean vector with components x and y, 

the unknown coordinates of the object, and is the covariance 

2 2 ... 
matrix with elements a.^ , a.^ and a • maximum likeli- 

IX ' ly ixy 

hood estimates x and y of the unknown coordinates x and y are 
the solutions of the two simultaneous linear equations deter- 
mined by 



9 ( log L) 
Tx 



x=x 

/N 

y=Y 



0 and 



9 (log L) 
9y 



/S 

x=x 

y=y 



0 



( 
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The equations can be written as 



A X + B y = D 

B X + C y = E 



and their solution as 



X = 



CD - BE 
AC - B^ 



y = 



AE - BD 
AC - B^ 



where A = Za., B = Zb., C = Zc., D = Z(a.x. + b.y.), 

1 1 1 11 i-'i' 

E = Z(b.x. + c.y.), a. = o.-' /d., b. = - o.^-v/d., c. = o-^'^/d. 

1 1 i-' 1 1 ly 1 1 ixy 11 IX ' 1 



2 2 2 

d. = 0 .^ o.^ - o-''^ and all of the sums are for i from 1 to n 
j. XX xxy 



since they are linear combinations of the estimates x^ and 

y^ , the estimates x and y are determined by a bivariate normal 

distribution. Consequently, all that is required to determine 

this distribution is its mean vector with components and 

2 2 

and its covariance matrix with elements and The 

X ' y xy 

mean vector is determined by 



y- = E{(CD - BE)/ (AC - B^)} = x 

and 

y- = E{(AE - BD)/(AC - B^)} = y 



And, the covariance matrix is determined by 

= E { (CD - BE) - E (CD - BE) }^/(AC - B^) ^ 

X 

= {C^ (F+I + 2L) - 2CB (H+K+M) + B^ (G+J+2N) }/ ( AC - B^)^ 
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= E{(AE - BD) - E (AE - BD) }^/(AC - B^)^ 

= {A^(G+J+2N) - 2AB(H+K+M) + B^ (F+I+2L) }/ (AC - B^ ) ^ 



2 > 2 



= E{ [ (CD-BE) -E (CD-BE) ] ♦ [ (AE-BD) -£ (AE-BD) ] }/ (AC - B ) 



= { (AC+B^) (H+K+M)-CB(F+I+2L)-BA(C+J+2N) }/(AC - B^)^ 
where F = , C = H = Za.b.a.^^ , I = Zb.^a.^'^, 

1 IX 1 IX 1 1 IX ' 1 ly 

2 2 2 2 
J = Zc . o.^ , K = Zb.c.o.-' , L = Za.b.o.^^, M = Z(a-c.+b. )o.^^ 

1 ly 1 1 ly ' 11 ixy i i i ixy 



and N = Zb.c.o.^-' where all of the sums are for i = 1 to n. 
1 1 ixy 



By using arguments given in Appendix 1 of this report, one 

can show that the axes of the elliptical confidence regions 

associated with x and y are coincident with an x ' y ' -coordinate 

system where the transformation from the xy-coordinate system 

to this system is the coordinate axes rotation through the 

2 2 

angle y defined by tan 2y = 2o"'^/(a^ - O'' ) . For a confidence 

p, from Reference 3, the minimum area confidence region is an 

2 

ellipse with semi-axes k and k > and area irk o^iC^i 

X y / X y 

where k = [-2 In(l-p)] ^ , > 



X * 


2 2 
= 

X 


2 

cos Y 


- 2o^ ^ 
xy 


COS 


Y 


sin 


Y 


+ 


2 

y 


. 2 
sin 


Y , 


< >1 
D 


2 2 
= 

X 


. 2 
sin Y 


+ 2o^ ^ 
xy 


cos 


Y 


sin 


Y 




2 

Y 


2 

cos 


Y 



/N 

and the center of the ellipse is at the point (x,y). In this 
coordinate system, o^iyi = 0 . 
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The above equations can be used to specify a composite posi 

tion estimate in terms of the location, orientation and size of 

an elliptical confidence region which is generally the form in 

which position estimates of the kind that are being considered 

2 2 

here are specified. But, since values of and 0^.9 

X y xy 

are required for each of the n estimates that is being combined 
a way is needed for determining these values given the orienta- 
tion and size of an elliptical confidence region. A procedure 
to do this when the orientation is given in terms of the direc- 
tion 6 of the major axis and the size is given in terms of the 
lengths SMJ and SMI of the semi-major and semi-minor axes and 
the confidence p is described next. 

By using an xy-coordinate system in which the positive 
y-axis direction is north and the positive x-axis direction is 
east and with the convention 0° < 6 < 180 °, the dependence of 
the value of the rotation angle Y and of the order relation 
between o^, and o^, on the value of the major axis direction 6 
is indicated by the following table; 



0° 


< 


6 


< 


45°: 


Y = 6 




and 


0 ^ , 

y' 


> 


Q 
X > 


45° 


< 


6 


< 


135° ; 


Y = 6 - 


0 

0 

(T\ 


and 


0 ^ , 


> 


0 ^ , 
















X * 




y' 


135° 


< 


6 


< 


180° : 


Y = 6 - 


180° 


and 


"y- 


> 


°x' 



With an order relation and a value for p, values for 0^, and 
o^, can be determined with values for SMJ and SMI. With values 

y 

2 2 

for O'', ,0^, and y, values for and 0-'-' can be found 

X ' ' y ' ' ' X ' y xy 

from the following equations: 
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2 2 2 2 2 

^X' ^ °y' 

2 2 2 2 2 

0 ^ = 0^, sin Y + o"'i cos y 

y X ' y 

and 

2 2 

0^^ = (o^, - 0^ ) sin Y cos y 

xy y ' X ' ' 

which can be obtained by inverting the equations above for o-' , 

, 2 
and 0 ^ 

y 

As an example, suppose the data in the following table 
represent three independent position estimates: 







X y 


6 


SMJ 


SMI 


k 


1st 




-3.7 18.1 


59° 


36 


20 


2 


2nd 




11.8 8.4 


105° 


37 


11 


2 


3rd 




0 0 


146° 


45 


23 


2 


Here 


, distances are in nautical miles and 


p = .86 in each 


case 


For 


this 


example with values in square nautical miles: 






2 

0, ^ 

lx 


= 264.58, 


= 159.42 


and 0 ^ 

Ixy 


= 98.89 






2 

Co 


= 321.35, 0.-^ 


= 51.15 


and 


= -78 






2y 


2y 




2xy 








2 

0,^ 


= 249.20, 0^-^ 


= 389.30 


and 


= -173.4 






3x 


3y 




3xy 






These values give the following composite 


estimate : 





X = -2.46 nautical miles and y = 12.26 nautical miles. 
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6.21 nautical 



For this case, = 10.28 nautical miles, = 

X y 

miles and = -55.1 square nautical miles. And, for k = 2: 

SMJ = 20.56 nautical miles, SMI = 12.42 nautical miles and 
6 = 119°. The composite confidence region and its three com- 
ponent confidence regions are shown in Figure 5. 

As a second example, suppose each position estimate is deter- 
mined by a circular normal distribution. Then a • = o.^ = o. 

IX ly 1 

and = 0 for i = 1 to n. In this case, the composite 

estimate is: 

X = (Zx./o.)/(Zl/o^) , y = (Iy^/a.)/(Il/a.) , 

= n/(El/a^)^ , = n/(El/o^)^ and = 0. 

In this example, since x and y are determined by a circular 
normal distribution, the minimum area confidence regions are 
circles and orientation is not an issue. 



N 
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Y 



North 




Figure 5. The ellipses define the confidence regions of 
the first example. The composite confidence region is in 
black. The position estimates are at the center of the 
ellipses. The numbers indicate the order of the estimate 
in the table on Page 17. 
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APPENDIX 1: A Position Estimation Procedure 



Consider a target whose position is unknown and a set of 
reference stations whose positions are known. Assume conditions 
are such that observed target bearings from the stations can be 
considered to be values of independent normal random variables 
whose means are equal to the true target bearings and unknown but 
whose standard deviations are known. The position estimation 
procedure that is described here is a maximum likelihood estimation 
procedure that is based on this assumption and the assumption that 
the conditions are such that the stations and the target can be 
considered to be located on a plane tangent to the earth's 
surface at a point in their neighborhood. 

With the above assumptions, the likelihood of observed 
bearings 0^^, 02 / •••» 0j^ from stations labeled 1, 2, ..., n is: 



L(0^, 02, ..., 



n exp 

1 /Tn e . 

1 



1 

2 



n 

Z 

1 



(0 



A ^ 2 / 2 

. ) /e . 
1 1 



where <t> 2 > •••/ are the unknown station true bearings and 

e,, e~, ..., e are the known stations standard deviations. 

1' 2 n 

To a first order approximation, the set of bearing estimates 

/\ /V /N. 

(Jf, (p^, ..., <I)^ that are determined by the procedure make 

L(0^, 02, ..., 0^) a maximum subject to the constraint that the 

bearing lines determined by a set of bearing estimates must all 
pass through a common point. The common point is the estimate 
of the target's position. 
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Figure 




6. The coordinate geometry. The coordinates of the initial 
estimate are (x*,y*). In the development, the reference 
point is at the initial estimate. 
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In order to impose the constraint on the estimate through 

consider the quantities u^ = = ^i ^^i ~ ^ 

where i represents any station number from 1 to n. 

Here, is the target's true bearing, (i^ is the bearing of an 

initial estimate of the target's position, is an observed target 

bearing and r^ is the range of the initial estimate from the 

station. It is also the radius of a circle that is centered on 

the station and passes through the initial estimate as shown in 

Figure 6 and u^ , and w^ are related arc lengths on this circle 

with u. = w. - V. . In this relation, w. is known and v. can be 
111 1 1 

expressed in terms of the unknown coordinates of the target with 
an approximation that does not involve i|)^. To do this, consider 
a rectangular coordinate system whose origin is at the position 
of the initial estimate and whose axes are oriented like those 
shown in Figure 6. To first order in this system, v^ = x cos - 
y sin where x and y are the unknown target coordinates and 
0 .-(}). = (0. - B.) - (x cos B- - y sin 3.)/r. . The use of this 
relation implies that the bearing line determined by B^, the 
bearing of the initial estimate, is approximately parallel to the 
true bearing line determined by The use of this relation for 

all stations imposes the constraint on the maximum likelihood bearing 
estimates by replacing what would otherwise have been estimates 
of n independent bearings 4']^ through (})^ by estimates of two independent 
quantities, the rectangular coordinates x and y. It also implies 
that the initial estimate's range from a station is approximately 
the target's range from a station, that is, that the initial 
estimate's position is relatively close to the target's position. 
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I 



since (0^ - 4>^) , the likelihood of observed bearings 



9,, 9„, .... 0 can be written as: 

1 2 n 



L (f 



’l' °2' * 






1 

n = — 

1 /2tT a. 

1 



1^2 2 
exp - Tr Z u./a. 

2 1 1 



where o. = r.e. and e. is the standard deviation of 0.. The 
1111 1 

maximum likelihood estimates for 6, through b are determined 

1 n 

by the estimates for x and y that make L(0j^, Q^, ..., 9^) a 

In this case, making L(9,, 9„, 9 ) a maximum is 

^ 12 n 

2 2 

equivalent to making 1 (u./o.) a minimum. So, to find the maximum 

1 ^ ^ 

likelihood estimates x and y, solve the following two equations 
for X and y : 



maximum. 



3 (tn L) 
3x 



= 0 



and 



3 ( tn L] 
dy 



= 0 



X = X 



X = X 



y = y 



y = y 



with the constraint given by u^ = w^ - x cos 3^^ + y sin 3^^ where 

w. = r. (9. - 3.) , the equations can be written as follows: 

1 11 1 ^ 



^ . 2 
Z (w. - X cos 3. + y sin 3.) (cos 3 . ) /o . = 0 
1 1 ■' 1 11 



and 



Z (w. - X cos 3. + y sin 3.) (sin B.)/o. = 0. 

2^ i 1 1 11 
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In terms of the following quantities: 



A = Z(cos^ , B = Z(sin 6^^ cos 6^)/o? , 

C = Z(sin^ B^)/o^ , D = L(w^ cos 6^)/a^ , 

E = Z (w^ sin , 

the equations become: 

Ax - By = D 

Bx - Cy = E 

The solutions are: 

(1) X = (BE - CD)/{B^ - AC) 
and 

(2) y = (AE - BD)/(B^ - AC) 

A confidence region can be constructed about an estimated 
position. In order to indicate how this can be done, a proba- 
bility region about the true position will be considered first. 

Both X and y are values of random variables. If a new set 
of bearings 0^^, © 2 / observed (for the same initial 

estimate and a fixed target), in general, a new pair of values 
X and y will be obtained. 
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If X and Y represent these random variables, from (1) and (2) , 



( 3 ) 



1^2 

X = — 5 Z (W./o,*) (B sin 3,- - C cos 3.) 

(b"^-AC) 1 ^ ^ ^ ^ 



(4) 



^ 1^2 

Y = — = Z (W./a.) (A sin 3- - B cos 3-) 

(B -AC) 1 ^ ^ 



where W. = r.(0. - 3-). 

1 11 1 

Since X and Y are a linear combination of the n normal random 
variables , W 2 , W^, or equivalently of the n normal random 

variables 0^^, © 2 , 0^, they have a joint normal distribution. 

Since E(W^) = ~ , if 3^ = <j)£ for i = 1, 2, n, that 

is, if the initial estimate of the target's position is at the 
target's position, E(Wj^) = 0 for i = 1, 2, ..., n. In this case 
E(X) = 0 and E(Y) = 0 and the joint normal distribution is centered 
on the object's position. To the degree of the approximations 
that have been made, this is also true if the initial estimate is 
not at the target's position. 

A region of minimum area for a given probability of contain- 
ment of an estimated position can be determined. The region is 
bounded by an ellipse which is centered on the object's position 
and whose axes lie along the axes of an x ' y ' -coordinate system 
that is obtained by rotating the xy-coordinate system that is 
centered on the object's position through an angle y. In this 

/s /s 

system, is 0, that is X' and Y' are independent normal ran- 

dom variables. The two coordinate systems are illustrated in 
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Figure 7. The coordinates of a point in the two systems 
related by 



are 



x' = X cos Y - y sin Y 

y' = X sin Y + y cos Y 

These relations imply: 

2 2 2 2 2 

(5) cos Y “ 2a cos y sin y + o sin y , 

XX xy y 

,^,22.2 22 

(6) sin Y + cos y sin y + cos y 

y' X y 

and 

2 2 2 2 

(7) ^ = (0^ - 0 ^) sin Y cos y + (cos y - sin y) 

x'y' X y xy 



where y, the angle of rotation of the coordinate axes, is posi- 
tive in the clockwise direction. And = 0 implies 



2 a^ ^ 

tan 2 y = ~ 2 — 2 

y X 

with the initial estimate of the target's position at the 

2 

target's position E(W^) = 0 and therefore Var(W^) = for 

i = 1, 2, n. In this case , from (3) and (4) 



2 

X 



2 2 
(B^-AC) 



2 

Z (l/o . ) (B sin 6 . 
1 



C cos 



6i) 



2 
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North 




I East 

Figure 7. An elliptical confidence region and the primed 
coordinate system in which the covariance 
is zero. The center of the ellipse and the 
origin of the coordinate systems are at the 
target's estimated position. The estimated bear- 
ing a and estimated range p are indicated for a 
reference point. The major axis direction is 6. 



28 



1 



2 

o 

/N 

y 



1 " 2 2 
— 5 y Z (1/a.) (A sin 6. - B cos 6.) 

(B^-AC) 1 ^ ^ ^ 



and 



12 . 

^ - — 5 9 ^ (l/a.)(B sin 6. - C cos 6.)(A sin 6. - B cos 3-)- 

xy (B^-AC)^ 1 ^ ^ ^ ^ ^ 



Using the definition for A, B and C, the above become 



( 8 ) 

(9) 



and 



„2 _ C 

° 2 ' 
X (AC-B ) 



^ 2 ' 
y (AC-B ) 



( 10 ) 



a 

/N /s 

xy 



(AC-B^) 



So tan 2Y = 2B/(A-C) for 6^^ = i = 1, 2, n. 

With the target's position known and, consequently, ())^ known 

2 2 

for i = l, 2, ..., n, the above expressions for and y 

f X y ' xy 

can be used, since the initial estimate of the target's position 

can be taken as the target's position. 

With values for and y, values for a^, and 

X y xy X ' y ' 

can be found by using equations (5) and (6). The probability 
that an estimated position will be within an ellipse of semiaxes 
ko'', and ka^ , which is centered on the target's position is 
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2 

1 -exp(-k /2) . This result can be found by integrating the bi- 
variate normal density over the ellipse. And the area of the 
2 

ellipse IS k 

t' x' y ' 

Given estimates x and y found by using Equations (1) and (2) , 
the ellipse with semi-axes and ko^, in a coordinate system 

/S /N 

that is centered on the point (x,y) and has been rotated through 

2 

an angle y is a 1 - exp(-k /2) confidence region. This follows, 
since, to the degree of the approximations involved, the bivar- 
iate normal distribution of X and Y is centered on the target's 

. 2 2 
position. The confidence ellipse is defined if and o^-' 

^ ^ X ' y xy 

can be found, that is if the elements of the covariance matrix 

can be found. To the degree of the approximations involved, 

this can be done as follows: First, assume the initial estimate 

of the target's position is at the target's position. Then, 

2 2 

values for ^ , “^xy determined by using Equations 

(8) , (9) and (10) These values can then be used to determine 

2 2 

°x' ' "^y* ^x'y' using Equations (5), (6) and (7). Now, 

with a value for k, a confidence region can be constructed. To 
the degree of the approximations involved, the shape of the con- 
fidence region is independent of both the target's position and 
of the initial estimate of the target's position. 

For the case where bearings are taken from the target on two 
or more stations, 0^ is the reciprocal of the bearing taken from 
the target. 

A discussion of the theory of bearings only position 
estimation procedures for situations similar to the one considered 
here is given in Reference 2. 
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The following equations are evaluated in the program to 
determine the coordinates x* and y* of the initial estimate: 



X* sin sin (a^-9^) ] sin 6^ 

- L p 2 sin (a 2 - 02 )] sin 6^ 



and 



y* sin = L Pj^ sin cos 02 

- L P 2 sin (a2~02)] cos 0^^ . 



Using the point determined by two lines of bearing as 
the initial estimate was suggested by a similar procedure 
described in Reference 3. 
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APPENDIX 2; 



Program Descriptions 
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PEST is a program that implements a position estimation 
procedure that is described in Appendix 1. With two or more 
observations on or from an object from or on two or more 
stations, the program will generate a position estimate for the 
object with an associated confidence region. After program 
initiation, the program user is presented with two options; (1) 
input bearings from stations on an object, (2) input bearings on 
stations from an object. To choose the first option, press T 
or t. To choose the second option, press S or s. Bearing 
information is entered in the following order: (1) the observed 

bearing on or from a station, (2) the bearing of the station 
from a reference location, (3) the range of the station from the 
reference location and (4) the bearing error. The reference 
location can be a station location. In this case, the station's 
bearing and range from the reference location are both taken to 
be zero. After two observed bearing observations and their 
associated station and bearing error information have been input, 
the user is given the option of (1) generating an estimate or 
(2) continuing to enter bearing observation data. A position 
estimate for an object is specified in terms of its range and 
bearing from the reference location. An associated elliptical 
confidence region which is centered on the estimate is specified 
in terms of the length and direction of its major axis, the 
length of the its minor axis and the probability (confidence) 
that it contains the object. The program user is given the 
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option of defining the confidence region in ternis of (1) size 
or (2) containment probability. For example, the size of an 
ellipse with a containment probability of .8647 is 2. In this 
case, the ellipse is referred to as a two sigma ellipse, since 
the semi-major and semi-minor axes of the ellipse are two standard 
deviation units in length. To choose the first option, press 
S or s. To choose the second option, press P or p. The estimate 
can be recalled by pressing E or e. Additional observations can 
be entered by pressing C or c. To quit, press Q or q. 

The position estimation procedure requires an initial 
estimate of the object's position. In the program, the initial 
estimate is at the intersection of the two bearing lines that 
correspond to the first two bearings that are input to the 
program. Because of this, the first two bearing inputs should 
be from or on the two stations that are estimated to provide the 
best initial estimate. Note, if bearing errors are large 
relative to the angular separation of the two stations as 
seen from the object, the bearing lines from the two stations 
may not intersect. If this is the case, the reciprocal bearing 
lines will intersect and a gross error in the final position 
estimate could occur. 

The procedure which is essentially equivalent to one 
described in Reference 2 is based on the following assumptions: 

(1) Bearings are taken on or from an object of unknown position 
from or on two or more stations of known position. (2) The 
distances involved are such that the object and the stations 
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can be considered to be located on a plane surface (a flat earth) . 
(3) The values of observed bearings on or from stations are 
determined by independent normal random variables each with a 
known standard deviation (the bearing error) and a mean equal 
to the true bearing (zero bias) . (If there is bias, it is known 
and removed.) 
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COMP is a program that implements a procedure for combining 
a group of position estimates for an object. In order to use the 
procedure, the estimates must have associated elliptical 
confidence regions. When used, the procedure combines the 
estimates and their associated confidence regions into a single 
composite position estimate and associated elliptical confidence 
region. After program initiation, the program user is prompted 
to input the number of elliptical areas (confidence regions) to 
be combined. The user is then presented with two options for 
specifying the areas: (1) by containment probability or (2) by 

size. For the first option, the areas are specified by the 
probability (confidence) that they contain the object. This 
option is chosen by pressing P or p. For the second option, the 
areas are specified in terms of ellipse size. This option is 
chosen by pressing K or k. The size of an ellipse with a 
containment probability of .8647 is 2 . In this case, the ellipse 
is referred to as a two sigma ellipse, since the semi-major and 
semi-minor axes of the ellipse are two standard deviation units 
in length. The orientation of an ellipse is specified in terms 
of the direction (the angle delta in the program) of its major 
axis. The position estimates, the centers of the elliptical 
areas, are specified in terms of latitude and longitude. The 
user is given the option of determining a confidence region in 
terms of (1) its containment probability (confidence) or (2) its 
size. To choose the first option, press P or p. To choose the 
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second option, press K or k. To recall the position estimate, 
press E or e. To quit, press Q or q. 

The procedure which is described in detail in Section III is 
based on the following assumptions: (D The individual estimates 

are determined by independent normally distributed random 
vectors with known covariance matrices and a common but unknown 
mean vector whose components are the object's rectangular 
coordinates. (2) The object is located in the plane of the 
coordinate axes, a plane tangent to a spherical earth. In 
the program, the coordinates of the first entry are the coordi- 
nates of the point of tangency. Because of other uncertainties, 
this flat search assumption should not introduce significant 
estimation errors for the distance scale for which the program 
is intended. In a sense, the procedure is a generalization of 
the procedure described in Reference 4 which is limited to 
circular confidence regions. 
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APPENDIX 3: 



Program Listings 
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PEST 



PRINT "S=bearings on the 



DR A:^="t" THEN JJ=0: 
OR A:$="s" THEN JJ=1: 
"observed bearing";P! 



GOTO 60 

GOTO 60 ELSE GOTO 25 
P=P*ATN(l>/45: IF JJ=1 



THEN 



FOR M=0 
P=A (M) : 



TO 1 

Q=A<M+2) : 



"station 



IF 1=1 THEN 



Z=SIN(A( 1) -A(0) 



A(0) ) 



R=A(M+4>; 0=A(M+6); GOSUB 390; NEXT M: GOTO 



10 CLS: SCREEN 0,0: DIM A (7) 

20 PRINT "T=bearings on the target. "i 
stations. " 

25 A#= INKEY# 

30 IF A#="T" 

40 IF A#="S" 

60 CLS: INPUT 
P=P+180 

70 INPUT "station bearing";Q: Q=Q*ATN ( 1 > /45; INPUT 
range";R; INPUT "bearing error";0: 0=0*ATN ( 1 > /45 
80 IF i=2 THEN GOTO 140 

90 1=1+1: A(I-1>=P: A(I+1)=Q: A(I+3)=R: A(I+5)=0: 

GOTO 60 

100 X=A(4)*SIN(A(2)-A(0) > : Y=A (5> *SIN (A (3) -A ( 1 > ) : 

) ; IF Z=0 THEN GOTO 330 

110 U=(X+SIN( A(1>>-Y+SIN (A(0>>)/Z; V=(X*COS ( A ( 1 ) ) -Y+COS ( 

) /Z 
120 
130 
150 

140 GOSUB 390 

150 CLS: PRINT; PRINT "E=Est C=Cont" 

160 A#= INKEY#: IF A#="E" OR A#="e" THEN GOTO 180 

IF A#="C" OR A#="c" THEN GOTO 60 ELSE GOTO 160 
CLS: F=(B*B-A*C): IF F=0 THEN GOTO 380 

X1=U+(B*E-C»D) /F; Y1=V+ ( A+E-B+D) /F; GOSUB 440: K=R1: J=B1 
T=SGN (B)*ATN( 1 ) : IF A=C GOTO 220 

T=. 5*ATN (2*B/ (A-C) ) 

G=<C*C0S(T)*C0S(T)-2*B*C0S(T)*SIN(T)+A*SIN(T>*SIN(T> ) /-F: 
G=SQR(G) 

230 H= (C*SIN (T) *SIN (T) +2*B*C0S (T) *SIN (T> +A-K-COS (T) *COS (T> > /-F: 

H=SQR(H>; IF H>=G GOTO 250 

240 Z=H; H=G: G=Z: T=T+2*ATN ( 1 ) 

250 CLS; PRINT: PRINT USING "\ \####. ##" ; "bearing=" , J*45/ 

ATN( 1 > 

255 I1=K: GOSUB 500: PRINT "range=";Il; PRINT 

260 PRINT "S=Size P=Prob E=Est C=Cont Q=quit" 

265 A#= INKEY# 



170 

180 

190 

200 

210 

220 



270 

280 

290 

300 

310 

320 

330 

340 

370 

350 

360 

370 

371 

372 
375 
X*Y; 



IF A#="S" 

IF A#="P" 

IF A#="E" 

IF A#="C" 

IF A#="Q" 

INPUT "size";S: 
0=l-EXP(-S+S/2> 
PRINT USING "\ 



OR A#="s" 
OR A#="p" 
A#="e" 
A#="c" 
A#="q" 



OR 

OR 

OR 



THEN 

THEN 

THEN 

THEN 

THEN 



CLS; 

CLS: 

CLS; 

CLS: 

END 



GOTO 

GOTO 

GOTO 

GOTO 

ELSE 



320 

350 

250 

60 

GOTO 265 
320 



IF S<=0 THEN GOTO 

\##. ####"; "probabil ity=" ,0; GOTO 



INPUT "probabil ity";0: IF 0>=1 OR 0<=0 THEN GOTO 350 

S=SQR(-2*L0G(l-0) > : PRINT USING "\ \####. ##" ; "si ze=" ,S 

X=S*G: I1=2*X: GOSUB 500:PRINT "major axis=";Il 

N=T*45/ATN(1) : IF N<0 THEN N=N+180 

PRINT USING "\ \####.##"; "direction=",N 

Y=S*H; Il=2+Y: GOSUB 500: PRINT "minor axis=";Il; I1=4*ATN(1>* 
GOSUB 500; PRINT "area=";Il: GOTO 260 



380 PRINT "no solution": END 
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90 X1=U-R*SIN(Q) : Y 1 =V-R*C0S ( Q) : GO 3UB 440 
00 W=P-B1: L=R1*0: IF L=0 THEN GOTO 380 

10 G=C0S(B1>/L: H=SIN(B1)/L; IF W>=4*ATN(1) THEN W=W-8*ATN ( 1 ) : 
OTO 43'> 

20 IF W<=-4*ATN(1> THEN W=W <-8*ATN ( 1 > 

30 W=W/0; A=G*G+A: B=G*H+B: C=H*H+C: D=W*G+D; E=W*H+E: RETURN 
40 R1=:3QR(X1*X1+Y1*Y1> ; IF R1=0 THEN B1=0: RETURN 
50 IF (^BS(X1/R1)=1 THEN Ml =SGN ( X 1 ) *ATN ( 1 ) *2 ELSE Ml =ATN ( X 1 /R1 / 
QR(1-X L*X1/R1/R1> > 

60 IF ABS(Y1/R1)=1 THEN BI=2*ATN ( 1 ) * (1-SGN ( Y1 ) ) ELSE B1=2*ATN(1) 
ATN(Y1/R1/SQR( 1-Y1*Y1/R1/R1 ) ) 

70 IF MKO THEN B1=8*ATN ( 1 > -B1 
80 RETURN 

00 11=100*11: I1=INT(I1): 11=11/100 

10 RETURN 
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COMP 



5 CLS: PI=4*ATN(1): INPUT "Number of Elliptical Areas"; AO: PRINT 

10 PRINT "Area Def i ni t i ons: PRINT 

11 PRINT "By Containment Probability, Press P.": PRINT 

12 PRINT "By Sigma Size, Press K. " 

15 E$=INKEYJ: 

16 IF E#="P" OR E^="p" OR ES="K" OR E:f="k" THEN GOTO 17 ELSE GOTO 
15 

17 CLS: PRINT "For the Latitude and Longitude entry format, press 
F. Otherwise, press C. " 

18 K^= INKEYS 

19 IF KS="C" OR KS="c" THEN GOTO 26 

20 IF KS="F" OR KS="f" THEN GOTO 23 ELSE GOTO 18 

23 CLS: PRINT "Latitude and Longitude are entered in degrees and 

minutes and tenths of minutes in the form: DDD-MM.MX where X is 

N,S,N or E. Leading zeros are optional, but the - must be 
included. Press C to continue." 

24 CS= INKEYS 

25 IF CS="C" OR CS="c" THEN GOTO 26 ELSE GOTO 24 

26 CLS: FOR 0=1 TO AO 

27 EE=0: PRINT 

28 INPUT "LAT";AS: GOSUB 410: IF EE=1 THEN GOTO 27 

29 Y=VA 

30 INPUT "LONG"; AS: GOSUB 410: IF EE=1 THEN GOTO 27 

31 X=VA 

32 IF 0>1 THEN GOTO 34 

33 Y0=Y: X0=X: C0=C0S ( Y0*PI / 180) 

34 Y=(Y-Y0)»60: X=(X0-X) 

35 IF X>270 THEN X=X-360 

36 IF X<-270 THEN X=X+360 

37 X=X*60*C0 

38 INPUT "Delta";S 

40 IF S>=180 OR S<0 THEN GOTO 35 

45 INPUT "Major Axis" ;P: INPUT "Minor Axis";Q: P=P/2: Q=Q/2 

50 IF ES="K" OR ES="k" THEN GOTO 60 
55 INPUT "P";P3: R=SDR (~2*L0G ( 1-P3) ) : GOTO 70 
60 INPUT "K";R 

70 IF S>=135 THEN S=S-180:T=P :P=Q: Q=T: GOTO 110 

80 IF S>=45 THEN S=S-90: GOTO 110 

90 IF S>=0 THEN T=P: P=Q: Q=T: GOTO 110 

110 P=P/R: Q=Q/R 

120 S=S*PI/180 

1 30 U=P*P#C0S (S)*COS(S)+Q»Q*SIN(S)*SIN(S) 

140 V=P*P*SIN (S) *SIN (S) +Q*Q*C0S (S) *C0S (S) 

150 W=(Q*Q-P*P)*SIN(S)*COS(S) 

160 Z=U*0-W»W: A=A+V/Z: B=B-W/Z; C=C+U/Z: D=D+V/Z*X-W/Z*Y 

170 E=E-W/Z*X>U/Z*Y: F=F+0/Z*0/ Z*U: G=G+W/Z*W/Z*U 

180 H=H-V/Z»W/Z»U: I=I+W/Z*W/Z*V; J=J+U/Z»U/Z»V 

190 K=K+W/Z*U/Z»0: L=L-0/Z*W/Z*W 

200 M=M+(V/Z*U/Z+W/Z*W/Z)*W: N=N-W/Z*U/Z*W 

210 NEXT 0 

215 CLS 

220 0=A*C-B*B: X= (C*D-B*E) /O: Y= ( A*E-B*D) /O 

230 R=(C*C*(F+I+2*L)-2*C*B*(H+K+M)+B*B*(G+J+2*N) ) / (0*0) 
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|>40 S= (A*A* (G+J-*-2*N) -2*A*B* (H+K+M) +B*B* (F+I+2*L) ) / (0*0) 

|>50 T= ( (A*C+B*B) * (H+K+M) -C*B* (F+I+2*L) -B*A* ((3+J + 2*N) ) / (0*0) 
i>55 Y=Y0+Y/60: VA=Y: GOSUB 470; GOSUB 500 
>57 PRINT "LAT = " ; 

>58 X=XO-X/ (60*00) ; VA=X: GOSUB 470: GOSUB 540 
l>59 PRINT "LONG = " ; G* 

'>70 A=SGN(T)*ATN(1) : IF S=R GOTO 280 

1>75 A=.5*ATN(2*T/ (S-R) ) 

|>B0 B=R*COS(A)*COS(A)-2*T*COS(A)*SIN(A) +S*SIN(A) *SIN(A) 

i>90 C=R*SIN(A) *SIN(A) +2*T*SIN (A) *COS (A) +S*COS(A) *COS(A) : A=A*180/ 

f'l 



500 

l>05 

>10 

j>20 

1525 

pO 

35 

;>40 

pi 

;>45 

>50 



GOTO 350 



>60 

f65 

j>70 

i>80 

>90 

)00 

no 

H5 

)20 

)23 

)24 

IE=1 

)25 

)30 

)35 

)40 

)50 

)60 

)65 



IF B>=C THEN D=A+90: GOTO 350 
E=B: B=C: C=E: IF A<0 THEN D=A+180: 

D=A; GOTO 350 

PRINT: INPUT "P";P3: K=SQR (-2*L0G ( 1-P3) ) 

PRINT: INPUT "K";K: P3= 1-EXP ( -K*K/2) : FI 

CLS: PRINT "Major Ax i s = ";2*K*SQR(R) 
PRINT "Minor Axis = ";2*K*SQR(S) 

PRINT "Delta = ";D; IF Fl=l THEN GOTO 345 
PRINT "P = ";P3: GOTO 350 
PRINT "K = ";K 
PRINT: PRINT "E = Est P 
A.f= INKEY# 

IF A#="E" 

IF A#="P" 

IF A#="K" 

IF A#="Q" 

CLS; PRINT 
CLS: END 
VA=VAL(A#) 



Fl=l; GOTO 330 



O 



= P K = K Q = Quit" 



OR A#="e" 


THEN 


GOTO 


390 




OR A#="p" 


THEN 


GOTO 


320 




OR A#="k" 


THEN 


GOTO 


325 




OR A#="q" 


THEN 


GOTO 


400 ELSE GOTO 


"LAT = "; 


L#: 


PRINT 


'•LONG = 


";G#; 


,4,1):: IF 


D#=" 


-" THEN GOTO 


435 


,3,1): IF 


D#=" 


-" THEN GOTO 


430 


,2,1): IF 


D#=" 


-" THEN GOTO 


425 


"Data entry error. 


restar t 


entry 


, 3 ) : GOTO 


440 








,4); GOTO 


440 








,5): GOTO 


440 









355 

GOTO 



550 



D#=MID#(A#, 
CLS: PRINT 
: RETURN 
M#=MID#(A# 
M#=MID#(A# 



OM=VAL(M#): VA=VA+VM/60 
R#=RIGHT#(M#, 1) 

IF R#="S" OR R#="s" OR R#="E" 
RETURN 



■for the area. 



OR R#="e" THEN VA=-VA 



170 QQ=SGN ( VA ) : VA=ABS ( VA ) 

)73 IF VA>180 THEN VA=(VA-360) ELSE GOTO 480 
)75 QQ=-QQ; VA=ABS(VA) 

mo DD=FIX (VA) : FF=VA-DD; MM=-FIX (6000*FF) / 100 
)90 RETURN 

00 IF QQ=1 THEN Q#="N": GOTO 520 
10 IF QQ=-1 THEN Q#="S" ELSE Q#=" " 

20 L#=STR# ( DD ) +STR# ( MM ) +Q# 

30 RETURN 

40 IF QQ=1 THEN Q#="W": GOTO 560 
50 IF QQ=-1 THEN Q#="E" ELSE Q#=" " 

60 G#=STR#(DD)+STR#(MM)+Q# 

70 RETURN 
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